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Abstract. Given a digraph G = (Va, Ag), a branching in G is a set of arcs B C Ag such 
that the underlying undirected graph spanned by B is acyclic and each node in G is entered 
(covered) by at most one arc from B. Tarjan developed efficient algorithms (based on the 
cycle contraction technique) for the following problem: given a digraph G with a weight 
function w: Ag —* M, find a branching B of the minimum weight w(B) := ^2 a£B u>(a) 
among all branchings with the maximum cardinality \B\. 

We generalize this notion as follows: for a digraph G and a matroid Mv on Vg, a matroid 
branching in G w.r.t. Mv is a branching in G such that the covered set of nodes is in- 
dependent w.r.t. Mv- The unweighted (cardinality) problem consists in finding a matroid 
branching B with |B| maximum. We show that the general cycle contraction approach is ap- 
plicable to this problem and leads to an efficient algorithm (provided that an oracle is given 
for testing independence in the matroids arising as the result of the contraction procedure). 
In the weighted version we are looking for a matroid branching B that minimizes w(B) 
(for a given weight function w : Ag — » K) among all matroid branchings of the maximum 
cardinality. We show that if Mv is a rainbow matroid (that is, nodes of G are marked with 
colors and it is forbidden to cover more than one node of any color), then there exists an 
0(min(n 2 , m log n)) method, matching the complexity of Tarjan 's algorithm (here n := \Vg\, 
m := \A G \). 

We also discuss a number of combinatorial tasks reducible to the weighted matroid branching 
problem. 
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1 Introduction 



For an arbitrary undirected graph G we write Vq (resp. Eg) to denote the set of nodes 
(resp. edges) of G. A similar notation is used for paths, cycles etc. If G is a directed graph 
we speak of arcs rather than edges and write Aq rather than Eq. For a set of nodes X 
denote the set of arcs entering (resp. leaving) X by 8 m (X) (resp. 5 out (X)). Denote the set 
of arcs having both endpoints in X by "f(X). 

A branching B in a digraph G is a set of arcs B C Aq such that the underlying 
undirected graph spanned by B is acyclic and each node in G is entered (covered) by at 
most one arc from B. Suppose that the arcs of G are endowed with real- valued weights 
w: Aq — > R. Then we get a problem as follows: 

(WB) Given G, w, find a branching B of the minimum weight w(B) among all branchings 
with the maximum cardinality \B\. 

(Hereinafter a real- valued function / : A — > R is assumed to be additively extended to 2 A 
by f(X):=Z x eAf^),XQA.) 

Chu, Liu [6j, and Edmonds [2] developed a polynomial algorithm for (WB) and also 
described the polytope V corresponding to the family of all branchings in B (that is, the 
convex hull of indicators of all branchings in G) . It turns out that V is given by the following 
set of inequalities: 

x: A G -> R + 

(1) x(5' m (v)) < 1 for each v G V G 

xll(X)) < \X\ - 1 for each X C V G 

Later, Tarjan [5] devised an 0(min(n 2 ,mlogn)) algorithm for (WB) (n := \Vq\, rn := 
\A G \). 

Recall [3J that a matroid is a pair (E,I), where E is a nonempty ground set (hereinafter 
assumed to be finite), X is a family of subsets of E such that: (i) G X; (ii) X £ I and 
Y C X imply Y G X; (hi) for each X,Y E X with \X\ < \Y\ there exists a G Y \ X 
such that X U {a} G X. The sets in X are called independent, others are called dependent. 
Inclusion-wise maximum independent sets are called bases. Moreover, for a given set X its 
inclusion-wise maximum independent subset is called a base of X. All bases of X are of 
the same cardinality, which is called the rank of X and is denoted by p(X). 

For a set B of arcs in G we denote by C(B) the set of nodes in G that are entered by 
an arc from B. For a digraph G and a matroid My = (T^j,Xy), a matroid branching B is 
a branching in G such that C(B) G Ty. We generalize (WB) as follows: 

(WMB) Given G, w, and My, find a matroid branching B of the minimal weight w(B) among 
all matroid branchings with the maximum cardinality \B\. 

We argue that (WBM) reduces to the weighted matroid intersection problem j4] and 
hence is polynomially solvable (if an oracle for testing independence in M.y is given). 
Indeed, M.y induces the matroid on Aq as follows: put Ma '■= (Aq,Ia), where X £ 1a 
iff each node in G is entered by at most one arc from X and C(X) G Ty. (The matroid 
axioms can easily be seen to hold for Ma-) 



We also consider the graphic matroid Mc '■= (Aq,Ic), where X G Z4 iff the underlying 
undirected graph spanned by X is acyclic. One can easily see that Xa n Iq is exactly the 
family of all matroid branchings in G w.r.t. My. Hence, polyhedral results [I] concerning 
the matroid intersection problem imply that the convex hull of indicators of all matroid 
branchings in G is given by 

x : A G -> R + 

(2) £ (x(S in (v)) : v G X) < p(X) for each X CV G 

x(l(X)) <\X\-1 for each X CV G 

where p denotes the rank function of My. This provides a natural generalization to (JTJ) . 

The rest of the paper is organized as follows. In Section [2] we consider the unweighted 
version of (WMB) (when w(a) = for all a), demonstrate the cycle contraction technique 
in action, and present an efficient algorithm based on it. Section [3] contains a discussion 
on so-called rainbow matroids and also describes a number of applications. An algorithm 
for (WMB) in case of a rainbow matroid is considered in Section |H Section [5] discusses 
implementation issues and, in particular, describes the 0(min(n 2 , mlogn)) algorithm that 
solves (WMB) for a rainbow matroid My. 

2 Unweighted Case 

Firstly, consider the following simplified version of (WMB): 
(CMB) Given G, w, and My, find a matroid branching B of the maximum cardinality \B\. 

To describe the algorithm for (CMB) we first introduce the notion of cycle contraction. 
For a given digraph G and a nonempty subset Q C V G construct the new digraph (denoted 
by G/Q) as follows: (i) all nodes in Q are removed; (ii) a new node (denoted by Q and called 
composite) is created; (iii) all arcs in j(Q) are removed; (iv) heads of all arcs in 5 m (Q) and 
tails of all arcs in 5 ont (Q) are changed to Q; (v) all other arcs remain unchanged. 

Since throughout the algorithm the nodes of digraphs are endowed with a matroid 
structure, one also needs a way to extend the notion of contraction to matroids. We consider 
an arbitrary matroid M = (E,T) and a nonempty subset Q C E. Put E' := E — Q U {Q} 
(where Q is a new composite element corresponding to the subset Q) and define the matroid 
(denoted by M/Q) on E 1 as follows: (i) a set X C E' — Q is defined to be independent 
w.r.t. M IQ iff there exists q G Q such that XUQ-{q} 6l; (ii) a set X C E' with Q G X 
is defined to be independent w.r.t. M/Q iff X — {Q} UQ Gl. 

Theorem 1. Let G be a digraph, My be a matroid on Vq, K be a cycle in G such that 
Z := Vk is independent w.r.t. My. Let B be a matroid branching w.r.t. My in G such 
that \BHA K \ = \Z\ - 1. Put G' := G/Z, M' v := My/Z. Then B' := B \ j(Z) is a 
matroid branching in G w.r.t. My. Moreover, if B is of the maximum cardinality, then 
so is B . 

Due to lack of space the proof of this statement is given in Appendix. 



Below we describe a generic version of the (CMB) algorithm. It has a certain amount 
of flexibility, which will be used later by the (WMB) algorithm. 

Algorithm. The algorithm for (CMB) consists of two phases: contraction and restora- 
tion. Let us describe the contraction phase first. The algorithm executes a series of steps 
as follows. At the i-th step it maintains the current graph G l , the current matroid My 
on V G i, and the current matroid branching B % in G l . Initially G° := G, M v := My, 
B° := 0. Consider the set U of nodes in G % that have at least one incoming arc in G l , put 
W := U - C{B i ). In case C^-B*) U {v} is dependent w.r.t. My for all v G W, it follows 
that B l covers a base of U . Therefore, B % is maximum; the algorithm stops. 

Otherwise, let v G W be a node such that C(B t ) U {v} is independent. Also, let a be 
an arc entering v. (In case of (CMB), these two choices are made arbitrary. Later on, we 
shall present an algorithm for (WMB) where a more careful selection of v and a will be 
necessary.) 

If B l U {a} spans no cycle, it forms a matroid branching in G % w.r.t. My. We proceed 
with the next step by putting 

G i+1 := G\ M$ l :=M v , B i+1 := B l U {a} . 

Otherwise, let K % be the cycle in B l U {a}. Denote by Z % the set of nodes of K l , put 

0*4-1 . = M^ 1 := M v /Z\ B l+1 := B i \ 7 (Z*) 

and proceed with the next step. 

The restoration phase gets the matroid branching B N in the final graph G N and re- 
verses the sequence of contractions (performed during the contraction phase) to obtain the 
matroid branching B° in the initial graph G° = G. Namely, suppose that cycle contraction 
took place on the z-th step; consider the composite node Z in G l . Two cases are possible: 

No arc from B l enters Z % . Since C(B % ) is independent w.r.t. M. l v , and ~Z? 
C(B l ) by definition of M v = My~ l /Z' l ~ l it follows that there exists a node w G Z 1 ^ 1 
in G 1 ' 1 such that C(B l ) U Z l ~ x — {w} is independent w.r.t. M l y X . Let a be the arc 
in K l ~ l that enters w. Now B % ~ 1 := B l U A K i-i — {a} is a branching in G l ~ x (hereinafter 
we identify arcs of G % with their images in G 1 ^ 1 ). Moreover, C(B t ~ 1 ) = C (B l )U Z 1 ^ 1 — {w} 
and hence B l ~ l is a matroid branching in G 1 ^ 1 . 

An arc b G B l enters Z l . Again, as C(B l ) is independent w.r.t. My and Z l G 

2 — 1 1 

C(B l ) it follows that C(B % ) U Z is independent w.r.t. My. Let a be the unique arc 
in K l ~ x that shares its head node with b. Now B l ~ x := B l U A K i-i — {a} is a branching 

in G 1 " 1 . Additionally, C(B i - 1 ) = (7(5^ - j^ 1 } U Z^ 1 and thus B*" 1 is a matroid 

branching in G 1 ^ 1 . 

The restoration phase terminates when reaching i = 0; B° is the required matroid 
branching in G° = G. One can easily see that B N is a matroid branching of the maximum 
cardinality in G N . Hence, from Theorem Q] and an inductive argument we get the following: 

Theorem 2. The above algorithm is correct, that is, B° is a matroid branching of the 
maximum cardinality. 



3 Rainbow Matroids 

It is clear that the above approach leads to a polynomial algorithm provided that an 
oracle is given that tests independence in matroids arising during the course of execution. 
Starting from here we only consider a certain special class of matroids admitting a simple 
independence test and closed under taking contractions (in the sense as we introduced it 
earlier) . 

For a nonempty set E the pair V(E) := (E,I), where 1 is the family consisting of 
the empty set and all singleton sets {e}, e G E, forms a partition matroid. For a pair of 
matroids Mi = {Ei,Xi), i = 1,2, with E\ f~l E 2 = 0, the direct sum of M\ and M 2 is the 
matroid Mi © M 2 = {E 1 U E 2 ,l), where 1 := {X U Y | X G X x , Y G 1 2 }. By a rainbow 
matroid we mean the direct sum of a number of partition matroids with disjoint ground 
sets. (The rationale for such a name is simple: the elements of the ground set of a rainbow 
matroid are assumed to be marked with colors; a subset is considered independent if it 
does not contain a pair of elements with the same color.) 

Let M = V(Ei) © ... © V{Ek) be a rainbow matroid with the ground set E := \J i Ei. 
For an arbitrary Z C E consider the matroid M' := M/Z. We claim that M' is also a 
rainbow matroid. More precisely, let / denote the set of indices i such that Ei n Z ^ 0. 
Then clearly 



That is, to contract an independent subset Z in a rainbow matroid we (i) "unite" colors 
of the elements from Z; (ii) remove Z from the ground set; (iii) add a composite element 
(that corresponds to Z and is denoted by Z) to the ground set. 

The (WMB) problem not only seems to be an appealing generalization to (WB), but 
also has a number of peculiar applications, which do not seem to be reducible to (WB). 
For example, let G be an undirected bipartite graph G with a bipartition X U Y of the 
nodeset Vq and a weight function w: Eq — > R on the edges. Additionally, suppose that the 
degrees of all nodes in Y do not exceed 2. Consider the problem of finding a matching M 
in G that has the minimum weight among all matchings with the maximum cardinality. 
We claim that this problem is reducible to (WMB) (for some rainbow matroid). Indeed, 
consider a node v G Y of degree 2; let e\ = {v,wi}, e 2 = {^,^2} be the edges incident 
to v. Construct the pair of nodes v\,v 2 and add the arcs a\ := (y\,v 2 ), a 2 := (v 2 ,vi) 
corresponding to e\, e 2 . Put w(ai) := w(ei), w(a 2 ) := w(e 2 ). For a node DGFof degree 1 
let e\ := {v,w\\ be the edge incident to v. Construct the pair of nodes v\,v 2 and the arc 
«i := (^1,^2) (corresponding to e±) with w(a±) := w(e±). The resulting graph is denoted 
by H. Finally, partition Vjj into color classes as follows. Each color corresponds to a 
node in X. A node v G Vh is marked with the color w G X if v has the inbound arc 
that corresponds to the edge of the form {v,w} G Eq, v G Y. The nodes v G Vh with 
5 m (v) = are assigned arbitrary colors. Clearly, there is a one-to-one weight preserving 
mapping between the set of matroid branchings in H and the set of matchings in G. 
As we shall show later, an optimum matroid branching in H can be found in 0(n log n) 



(3) 




time (where n := \Y\). This provides an improvement over the standard augmenting path 
approach, which leads to an Q{n 2 ) algorithm. 

The very same problem can also be restated as follows: given an undirected graph G, 
direct some edges in G such that: (i) each node of G is entered by at most one directed 
edge; (ii) the number of directed edges is as large as possible; (iii) ties are resolved by 
minimizing the total weight of all directed edges. Here each edge e € E G is assumed to be 
endowed with two reals indicating the weights for two possible ways of directing e. Using 
the approach as above, this problem can be solved in O(mlogn) time (where n := \V G \, 
m:=\E G \). 

4 Weighed Case 

We now proceed by presenting the algorithm for (WMB) in a graph G endowed with 
a weight function w: A G — > E and a rainbow matroid M.y = ©jf(Vi), Here V G is 
partitioned into the collection of sets {V\, . . . ,Vk} corresponding to different colors. Our 
goal is to find a matroid branching B minimizing w{B) among those with the maximum 
cardinality \B\. 

We say that a matroid branching B covers the i-th color if C{B) contains a node 
from V%. Firstly, we augment G to ensure that we are actually looking for a branching 
that covers all colors. To this aim, we add the auxiliary node s together with the auxiliary 
arcs (s,v) going to all other nodes. These new arcs are assigned large positive weights 
(e.g. max a w(a) + 1, where maximum is taken over all non-auxiliary arcs). Since no arc 
enters s, the color of s is irrelevant. 

Clearly, for the newly constructed graph there exists a matroid branching B covering 
all colors. Moreover, since the weights of the auxiliary arcs are large, B has the minimum 
weight in the augmented graph (among those covering all colors) iff the restriction of B 
to the initial graph yields a matroid branching of the minimum weight (among those with 
the maximum cardinality). Further on, we denote the resulting graph by G. 

Now the problem can be stated in terms of linear programming as follows: 

1 for each 1 < i < k 

for each X C V G 

-> min 

The correctness of this description can easily be derived from the standard polyhedral 
facts regarding the matroid intersection problem. Further, we shall present an algorithm 
that solves (WMB) and yields a 0, 1-solution to ([I]) (thus providing another proof for the 
correctness of the description). 

Consider the program dual to (|4|): 

vr: {l,...,k} -» E 
(5) £:2 v G^m + 

w' > 



x : A G -> R + 

E {x(5' m (v)) : v 6 Vi) = 

x( 7 (X)) < \X\ - 1 

Y\ (w(a)x(a) : a € A G ) - 



where 

w':=w-J2 <i) E ^ (V) + E * > 

i v&Vi X 

Here x A stands for the indicator of a set A. We shall never refer to the objective function 
of J5)) explicitly and hence we omit it for brevity. The weights w' are called reduced. 
Complementary slackness conditions for (SI) , dS]) are: 

(CS1) w'(a) = for each a € Ac such that x{a) > 

(CS2) x(j(X)) = \X\ - 1 for each X C V G such that > 

Algorithm. The algorithm for (WMB) is similar to that for (CMB) and also consists 
of the contraction and restoration phases. Contraction phase consists of steps. At the i-th 
step we maintain the current graph G 1 ', the current coloring of Vq%, the current matroid 
branching B l , and the current feasible solution (vr 1 ,^ 1 ) to j5]). Moreover, B % and (vr 4 ,^ 4 ) 
satisfy (CS1) (for tt := tt*, £ := f, a; := X^)- 

Initially G° := G, B := 0, £° = 0, and 7r° is chosen so as to satisfy w' > 0. At any 
step, the graph G l may be regarded as obtained from G° = G by contracting a number of 
disjoint inclusion-wise maximum sets Q%, . . . ,Q S C Vg (initially s = since no subset is 
contracted) . 

In case B % covers all the colors that have inbound arcs, we are done. Otherwise let k 
be an uncovered color that has at least one inbound arc. In contrast to (CMB) algorithm, 
in the weighted case two different kinds of steps are possible: primal and dual. If w' > 
for all arcs a entering nodes with the color k, then a dual adjustment is made by setting: 

n i+1 (k) := Tr^fc) +e 
(6) 7r i+1 (j) :=7T*(j) for all j ^ k 

(i i+1 (Qj)--=e{Qj) + £ for all j 

It is straightforward to verify that J6)) decreases w'(a) by e for each arc a entering a node 
with the color k while preserving all other reduced weights. The value of e is chosen as 
the largest positive number such that the reduced weights w' w.r.t. 7r l+1 and are 
nonnegative. Put G t+l := G 1 , B l+l := B % and proceed with the next step. 

A primal step (for the color k) is executed when there is an arc a such that: (i) a enters 
a node with the color k; (ii) w'(a) = 0. (In particular, by the choice of e as above we are 
guaranteed that each dual step is immediately followed by a primal one.) 

We try to add a to B l . If this does not lead to creation of a directed cycle, the i-th 
step is complete and we proceed with G' t+1 := G l and B l+1 := B l U {a}. Otherwise, let K 
be the cycle in B % U {a}; let Z be its nodeset. Contract K in G % by putting G l+1 := G l /Z 
(and updating colors according to j3j)). Also, put B l+1 := B % \ j(Z) and proceed with the 
next step. 

The restoration phase of the algorithm relies on bookkeeping from the contraction phase 
and essentially coincides with that in the algorithm from Section [2 

Theorem 3. The above algorithm is correct, that is, constructs a matroid branching that 
covers all colors and has the minimum weight. 



Proof. Since the algorithm for (WMB) essentially resembles the generic version for (CMB), 
it constructs a branching that covers all colors. It remains to show that this final branch- 
ing B has the minimum weight. To see this, we prove that B and the final functions 7r,£ 
satisfy (CS1), (CS2) (for x := x B )- By the standard linear programming arguments this 
implies the minimality of w(B). 

Let us call a subset Q C Vq contracted if it appears in the sequence Qi, ■ ■ ■ ,Q S on some 
step of the algorithm. First note that at the moment the algorithm contracts a cycle K one 
has w'(a) = for all arcs a of K . These arcs become hidden by contractions (fall into j(Q) 
for a contracted set Q) and the dual adjustments (J6j) never affect their reduced weights. 
Then, any current branching maintained by the algorithm in the current (contracted) 
graph obeys (CS1). Since the final branching B consists of arcs from the branching B N 
and certain arcs from contracted cycles, B also satisfies (CS1). 

As for (CS2): during the course of execution cycles are detected and contracted in the 
current graph. However, by taking preimage under contractions the nodesets of these cycles 
may be regarded as the subsets of Vq- Moreover, each of those nodesets is contracted. A 
simple inductive argument shows that \y{Q) C\ B\ = \Q\ — 1 holds for every contracted 
set Q. Also, it is clear that £(Q) > is only possible for a contracted set Q. Therefore, 
(CS2) is satisfied. 

5 Efficient Implementation 

Two efficient versions of the above algorithm can be presented achieving 0(m log n) and 
0(n 2 ) time bounds (here n := \Vg\, m := |^4g|)- We start with the first one. 

The graphs arising at the intermediate steps of execution are never stored explicitly. 
Instead the algorithm maintains the initial graph G and the collection of the contracted 
sets Qi, . . . , Q s . The Disjoint Set Union (DSU) data structure T> contr is used to store this 
collection (together with the singletons {v} corresponding to nodes v not covered by any 
of Qi). Here it is sufficient to use a version of DSU based on rank heuristic only [T] and 
thus achieving the O(logra) time bound for unions and root queries. The total number of 
unions performed by the algorithm is O(n), which totally results in the 0(n log n) term. 

The set of arcs of the current graph is also stored implicitly. That is, each arc of the 
current graph corresponds to a certain arc in the initial graph. However, not every arc from 
the initial graph survives the contractions. We call an arc from the initial graph dead if it 
is contained in one of j(Qi), 1 < i < s; other arcs are considered alive. To check if a given 
arc (u,v) is dead it is sufficient to compare the roots (with regard to V contr ) of u and v. 
Each such request is served in O(logn) time. The algorithm does not try to get rid of dead 
arcs as soon as they appear; instead the lazy cleanup strategy is used. The removal of a 
dead arc is postponed until this arc is discovered (see below). 

The current branching B is maintained by keeping, for each node v of the current graph, 
the arc from B entering v (if any). Recall that before adding an arc a the algorithm checks 
if BU{o} is acyclic. This operation is carried out as follows. The set of nodes of the current 
graph is partitioned into the equivalence classes: nodes u and v are considered equivalent 
iff they belong to the same directed tree of B. Clearly, this equivalence information can be 



maintained by another DSU instance T>t ree . Each time an attempt to insert an arc (u, v) 
is made, the algorithm checks (in O(logn) time) if u and v are equivalent. 

In case of the positive answer, B U {a} contains the cycle that can be extracted by 
following from v to the root of the corresponding tree. The time required to extract this 
cycle is proportional to its length; since the extracted cycle is immediately contracted, the 
sum of lengths of all extracted cycles is 0(n), which totally results in the 0{n log n) term. 

Otherwise (u and v are not equivalent) BU {a} is acyclic, so the algorithm proceeds by 
uniting the equivalence classes of u and v. Similar to previous, this yields the 0(n log n) 
term. 

Another DSU instance T> co \ is used to track the colors of nodes: each color c is repre- 
sented by the set of nodes of the current graph that are marked with the color c. Clearly, 
0(n log n) time is totally necessary to maintain this information. 

The dual variables ir are attached to the roots of T> co \. Variables £ are never maintained 
since they do not affect the reduced weights of the alive arcs. 

We now explain how reduced weights are maintained and how values e for the dual 
adjustments are calculated. The algorithm uses mergeable priority queues. A queue holds 
a set of {key, data) pairs. The keys are stored implicitly and may change during the course 
of execution. The key of a pair p is provided to the queue when p is inserted; each insertion 
costs O(logra) time. Later, the queue can report (in O(logn) time) the current value of the 
key in any pair. The queue can report (in (log n) time) a pair with the minimum key. Any 
pair that is currently stored in the queue may be deleted from it in O(logra) time. For an 
arbitrary real number 5, the queue can decrease by 5 the keys of all pairs contained in it. 
Finally, any two queues (71,(72 can be merged (in O(logn) time); this operation destroys qi, 
(72, and produces a new queue holding the union of sets that were held by qi, q%. A possible 
implementation for the described data structure is based on binomial heaps p]. 

Each color i is assigned the mergeable priority queue Q(i). This queue holds the set of 
pairs of the form (it/ (a), a), where a ranges over all arcs from 5 m (v ) and v is a node of the 
current graph that is marked with the i-th color. Additionally, due to the lazy nature of 
the cleanup strategy Q(i) may contain certain pairs corresponding to dead arcs. The keys 
of such pairs are of no meaning. To calculate the value of e for J6|) the algorithms makes 
the corresponding request to Q(k) (where k denotes the color that is chosen to be covered 
at this step) and extracts an arc a with the minimum value of key. Then it checks (in 
O(logn) time) if a is dead. If so, another arc is fetched. The procedure stops once an alive 
arc is discovered. Since each arc may be extracted at most once, the O(mlogn) bound for 
the total time spent for extractions follows. The adjustments (|6j) are executed by changing 
all keys in the corresponding queue by e; totally this takes 0(n log n) time. When a cycle 
is contracted and colors are united according to (|3j), the corresponding queues are merged. 
The sum of lengths of all contracted cycles is 0(n), which implies the 0(n log n) bound 
for all the merges. 

Summing up, we obtain the following: 

Theorem 4. (WMB) can be solved in 0{m log n) time. 

This algorithm can be slightly improved for the case of dense (m > n 2 /logn) graphs 
to obtain the 0{n 2 ) time bound. Firstly, we use simple (not involving any heuristics) 



implementations of DSUs T> contr , T> tree , and V co i . This way, a union operation is served 
in 0(n) time and a root query is answered in 0(1) time. In particular, the reduced weight 
of any arc can be found in 0(1) time. 

We represent the current graph by the matrix A, where A[u, v] holds the reference to 
the arc from a node u to a node v (if any). Then, to contract a set Z one needs 0(n \ Z\) 
time to update A, which totally results in 0(n 2 ) time for all contractions. 

Additionally, we maintain the matrix B. That is, for a node v and a color i, B[v,i] is 
an arc (if any) of the minimal reduced weight that leaves v and enters a node marked with 
the i-th color. To compute the value of e the algorithm scans the appropriate column of B 
and selects a minimum arc. Clearly, dual adjustments do change B. 

Contraction of a set of nodes Z is a two-step process. Firstly, the colors of nodes in Z are 
merged. To merge a pair of colors the algorithm scans the corresponding columns of B and 
chooses minima. This process takes 0(n\Z\) time, hence 0(n 2 ) time totally. The second 
step involves updates of B caused by merging nodes in Z. Clearly, only rows corresponding 
to nodes in Z are to be processed. These rows are replaced by a single row holding the 
minimum arcs going from the newly-created composite node Z. To this aim, we scan all 
pairs u £ Z, v ^ Z, look at the corresponding arcs A[u, v], and construct the Z-th row 
of B by taking minima. 

Theorem 5. (WMB) can be solved in 0(n 2 ) time. 
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A Proof of Theorem Q] 



Let Ma (resp. M' A ) be the matroid on Aq (resp. Aqi) corresponding to My (resp. M' v ). It 
follows that \C{B) n Z\ equals either \Z\ or \Z\ — 1. In the former case there exists a unique 
arc in B that enters Z. In the latter case no such arc exists. In both cases one can easily 
see that B' is a matroid branching in G' w.r.t. M' v . Moreover, C(B') = C{B) \ ZU {Z} 
in the former case and C(B') = C(B) \ Z in the latter. 

We construct the auxiliary bipartite digraph D with the nodeset Vd '■= Aq and the 
bipartition BUB, where B := Aq — B. For a pair of nodes x G B, y G B such that 
B — {x} U {y} spans no undirected cycle, we add the arc (y, x) to D. Similarly, for a pair 
of nodes x G B , y G B such that B — {x} U {y} is independent w.r.t. Ma, we add the arc 
(x, y) to D. An arc x G B is declared initial if -B U {x} is independent w.r.t. Ma and /ma/ 
if 5 U {x} spans no undirected cycle. A path in D going from an initial arc to a final one 
is called augmenting. It is known |1] that \B\ is maximum iff there is no augmenting path 
in D, 

Suppose that \B\ is not maximum and choose an augmenting path 

P = (ai, ...,a n ) 

of the shortest length in D, ai G Aq- Construct the digraph D' for G' and the matroid M' A 
(similar to D and Ma)- We argue that P can be transformed into an augmenting path P' 
in D' , which would imply a contradiction, as required. 

Suppose P contains an arc (x,y), x G B \ ~f(Z), y G B \ j(Z), By definition Bq := 
B - {x} U {y} is independent w.r.t. Ma- Since x <£ j(Z), then \C(B ) n Z\ > \Z\ - 1. 
Therefore, under the reduction from My to M' v the set C(Bq) generates an independent 
(w.r.t. My) set in Vc, which is equal to C(B' — {x} U {y}). So the arc (x,y) is also present 
in D'. _ 

Similarly, let P contain an arc (y,x), x G B \ "y(Z), y G B \ 'j(Z). By definition, y is 
contained in the fundamental cycle of B U {x} in G. This also remains valid for G' and 
B' U {x}, hence (y, x) is an arc of Z)'. 

There exists a unique node r G Z that is not covered by B n 7(Z). Let T denote the 
subtree of B rooted at r. Clearly, Z QVr and hence 7(Z) C 7(Vr). 

Let denote the smallest index such that G 7(Vy). (In case no such arc exists, none 
of a, is contained in j(Z) so it remains to check that the arc a\ (resp. a n ) remains initial 
(resp. final) in D' , see below.) Suppose cik G B, then a^-i G -B. Moreover, is contained 
in the fundamental cycle of B U {ak-i}, hence a^-i G j(Vt) — a contradiction with the 
minimality of k. Consequently, at G B. Suppose the head of is not r. Then, is not an 
initial arc since all nodes in Vr except for, possibly, r are covered by B. Now a^-i is the arc 
from B that shares its head node with at- Clearly at-i G ^(Vt), which again contradicts 
the choice of k. So we get that a& enters r. 

Let I denote the largest index such that ai G 'y(Z). The arc ai cannot be final, therefore 
I < n. One can easily see that a/ + i G B. Indeed, suppose the contrary. By definition of 
D', ai + i is contained in the fundamental cycle of B U {a/}. But a/ G j(Z), which implies 
ai + i G 'y(Z). This is, however, a contradiction with the way we choose I. 

Case splitting completes the proof as follows. 



Case 1: \C{B)nZ\ = \Z\ - 1. 

Subcase 1.1: /c = 1. We claim that ai + \ is an initial arc in D', hence 

P' ■= (ai+i, ■■■ ,a n ) 

is the required augmenting path for B' in D' . Indeed, B U {fi/+i} is dependent and B — 
{a/} U {a/ + i}, B U {a&} are independent (w.r.t. .Ma)- Therefore, by a simple reasoning 
one can show that Bq := B — {ai} U {a/^az+i} is an independent w.r.t. Ma set with 
\C{Bq) n Z\ > \Z\ — 1. Under the reduction from My to My, C(Bq) gives rise to an 
independent (w.r.t. My) set coinciding with C(B' U {a/ + i}). Therefore, the arc az+i is 
initial in D'. 

Subcase 1.2: k > 1. We claim that (afc_i,a/ + i) is an arc of D'. hence, 

P := (ai, . . . , ctfc-i, ai + i, . . . , a n ) 

is the required augmenting path for B' in D' . To see this, put Bq := B — {a^-i,^} U 
{ak,ai + i}. Since P has the shortest length, Bo is independent w.r.t. M.a (see [1]). One 
has \C(Bq) H Z\ > \Z\ — 1. Under the reduction from My to M' v , C(Bq) generates an 
independent w.r.t. Ai'y set, which coincides with C(B' — {ak-i} U {ai + i}). From this, the 
claim follows. 

Case 2: |C(i?)nZ'| = \Z\. We argue that similar to Subcase 1.2, (afc_i,a^ + i) is an 
arc of D' . Put Bq := B — {a/} U {a/+i} thus forming an independent w.r.t. Ma set. 
Clearly, \C{Bq) fl Z\ > \Z\ — 1. Under the reduction from My to My, C(Bq) generates 
an independent w.r.t. My set, which coincides with C(B' — {ofc_i} U {a/ + i}), as claimed. 

Finally, we deal with the arcs a\ and a n . The head node of a\ is not covered by B. If 
a\ £ 7(^) then |C(i?) PI Z\ = \Z\ — 1 and /c = 1. This is only possible in Subcase 1.1 but 
then P' starts with a/ + i, which exists in D' by the choice of I. Otherwise a\ £ 'y(Z) and 
the arc a\ remains existent in D' . Moreover, B U {a\\ is independent w.r.t. Ma, which 
implies that B' U {a±} is independent w.r.t. M' A , thus a\ is an initial arc in D'. 

Consider the final arc a n . The definition implies that B U {a n } does not span an 
undirected cycle in G. Hence, a n £ j(Z). Moreover, B' U {a n } still does not span an 
undirected cycle in G' , so a n is a final arc in D' . 

The proof of Theorem Q] is complete. 



